Anomalous Properties in the Normal and Superconducting States of LaRu3Si2 
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Superconductivity in LaRu3Si2 with the honeycomb structure of Ru atoms has been investigated. 
It is found that the normal state specific heat C/T exhibits a deviation from the Debye model 
down to the lowest temperature. A relation C/T = 7n + PT^ — ATlnT which concerns the electron 
correlations can fit the data very well. The suppression to the superconductivity by the magnetic field 
is not the mean- field like, which is associated well with the observation of strong superconducting 
fluctuations. The field dependence of the induced quasiparticle density of states measured by the 
low temperature specific heat shows a non-linear feature, indicating the significant contributions 
given by the delocalized quasiparticles. 

PACS numbers: 74.70.Dd, 74.25. Bt, 74.25.F-,74.10.-|-v 



Superconductivity arising from non-phonon mediated 
pairing, such as through exchanging the magnetic spin 
fluctuations, has renewed interests in condensed matter 
physics. The superconducting (SCj mechanism of the 
cupratesT] and the iron pnictidcs[2], although not yet 
settled completely, should have a close relationship with 
the electron correlations. 0-0] A similar assessment may 
extend to many others, like heavy FermionQ and or- 
ganic materials^]- In this regard, the systems i?r3Si2 
or RT^^2 {R stands for the rare earth elements, like La, 
Ce, Y, etc., T for the transition metals, like Ru, Co and 
Ni, etc.) provide an interesting platform, since a variety 
of combinations of chemical compositions allow the sys- 
tem to be tuned between superconducting and magnetic, 
and sometimes both phases coexist. [1, 0] Among these 
samples, the LaRu3Si2 has a SC transition temperature 
as high as 7.8 K.[10| The material of LaRu3Si2 contains 
layers of Ru with the honeycomb structure sandwiched 
by the layers of La and Si, forming a Pd^/ia or P63 
space group. Preliminary experiment found that the SC 
transition temperature drops only 1.4 K with the sub- 
stitution of 16 % La by Tm (possessing a magnetic mo- 
ment of about suggesting that the superconductiv- 
ity is robust against the local paramagnetic moment [8|. 
By doping the La with Gd, a coexistence of supercon- 
ductivity and the spin glass state ll| was observed. In 
CeRu3Si2, the SC transition temperature drops to about 
1 K and a valence fluctuation model was proposed for 
the pairing [1^. Since the Ru atom locates just below 
the Fe in the periodic table, a key player in the iron 
pnictide superconductors, therefore it is very curious to 
know whether the superconductivity here is induced by 
the electron-phonon coupling, or by other novel mecha- 
nism, such as the electron correlations. In this paper we 
report the results of transport and specific heat on sam- 
ples of LaRu3Si2. Our results reveal some novelties in 
both the SC and normal states of LaRu3Si2. 




FIG. 1: (color online) X-ray diffraction patterns of the sample 
LaRu3Si2. All main diffraction peaks can be indexed well 
by a hexagonal structure with a = 5.68 A and c = 7.13 A 
with Ru as the impurity phase. For some peaks the difference 
between the data and the fitting is a bit large because some of 
the grains of the polycrystalline sample are slightly oriented. 
The ratio between LaRu3Si2 and Ru is found to be 85:15. 
The inset gives a sketch of the structure. One unit cell is 
highlighted by the rhomboic block. 



The sarnples were fabricated by the arc melting 
method, d, [l^ Ell The starting materials La metal pieces 
(99%, Alfa Aesar), Ru powder (99.9%) and Si powder 
(99.99%) were weighed and mixed well, and presseed into 
a pellet in a glover box filled with Ar atmosphere (wa- 
ter and the oxygen compositions were below 0.5 PPM). 
In order to avoid the formation of the LaRu2Si2 phase, 
we intentionally let a small amount of extra Ru with the 
nominal compositions as LaRu3-|-a;Si2. Three round of 
welding with the alternative upper and bottom on the 
pellet was taken in order to achieve the uniformity. Af- 
ter these refined processes, the resultant sample contains 
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FIG. 2: (color online) (a) Temperature dependence of the DC 
magnetization measured in the ZFC mode and the FC mode 
at a magnetic field of 20 Oe. (b) The MHLs measured with a 
field sweeping rate of 50 Oe/s at different temperatures. At 
9 K, the MHL shows a rough linear paramagnetic behavior. 
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FIG. 3: (color online) Temperature dependence of resistivity 
at zero magnetic field. The inset shows the resistivity at dif- 
ferent magnetic fields: 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 
4.0 T. (b) Temperature dependence of the critical magnetic 
field with three different criterions: Hc2 (squares, 95%pn), Hc2 
(circles, 50%pn) and the irreversibility line Hirr (up-triangles, 
0.1%pn). There is a large area between the Hc2 (95%pn) and 
Hirr (0.1%p,i), which is probably induced by the strong SC 
fiuctuation. 



mainly the phase of LaRu3Si2 and small amount of Ru 
remains as the impurity phase. In Fig.l we plot the x-ray 
diffraction patterns (XRD) on one typical sample and the 
Rietveld fitting using the GSAS programfl^. It is clear 
that the main diffraction peaks can be indexed well by 
a hexagonal structure with a = 5.68 A and c = 7.13 
A. Some weak peaks arising from the impurity phase Ru 
can also be seen. A detailed fitting to the structural data 
find that the ratio between LaRu3Si2 and Ru is around 
85:15 for this typical sample. The sample preparation 
and the quality characterized by the SC transitions can 
be repeated quite well. It is found that, some of the 
LaRu2Si2 phase with a tetragonal structure can be found 
if the starting material has the nominal composition of 
LaRu3Si2. In this case, the XRD data exhibit clearly 
two set of structures and can be easily indexed by the 
GSAS program. For the present sample, the absolute 
difference between the experimental data and the fitting 
curve can be observed for some peaks because part of 
the grains in the sample are slightly aligned. The resis- 
tivity was measured with a Quantum Design instrument 
PPMS-16T with a standard four-probe technique, while 
the magnetization was detected by the Quantum Design 
instrument SQUID- VSM with a resolution of about 5 x 
10^^ emu. 



In Fig. 2(a) we present the temperature dependence 
of magnetization measured in the zero- field-cooling mode 
(ZFC) and the field-cooling mode (FC). By considering 
the demagnetization factor on the ZFC data, the Meiss- 
ner screening is estimated to be almost 100 %. This 
indicates that the SC connections between the grains of 
LaRu3Si2 are very good, although we have slight sec- 
ondary phase of Ru. The onset Tc determined from the 
magnetization is around 7.8 K. The majority of the SC 
transition occurs at about 6.6 K under a magnetic field 
of 20 Oe. This difference is not induced by the inhomo- 
geneity of the sample, it may be induced by the relatively 
strong SC fluctuations (see below). Fig. 2(b) shows the 
magnetization hysteresis loops (MHL) measured at dif- 
ferent temperatures. The symmetric and clear opening 
of the MHLs indicate that it is a type-II superconduc- 
tor. A roughly linear MHL was observed at 9 K, just 
above Tc, indicating that the normal state has no long 
range ferromagnetic order. We didn't observe a mag- 
netization enhancement near Tc, which was reported in 
Tm and Gd doped samples in early publications [sl. \v\. 
Fig. 3(a) shows the resistive transitions at zero field (main 
panel) and different magnetic fields under 4 Tesla (inset). 
The onset resistive transition temperature is at 7.9 K 
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FIG. 4: (color online) (a) The raw data of specific heat coef- 
ficient C/T vs. T^, at different magnetic fields ranging from 
to 5 T. The normal state data (at 5 T) shows a non-linear 
feature down to the lowest temperature here, indicating a de- 
viation from the Debye model, as shown by the red dashed 
line. The solid line represents the fit to the formula including 
the electron correlations (see text), (b) The electronic specific 
heat coefficient obtained by subtracting normal state value 
Cjv/T (data at 5 T) from the total. The solid lines are the 
theoretical fitting curves based on the BCS model, (c) The 
entropy difference (squares) between the SC state and the nor- 
mal state, derived from Ss - Sn = Joi^s/T' - Cn /T')dT' 
at zero field and 5 T. The condensation energy is calculated 
by Ec = Jo iSs-SN)dT'. (d) The magnetic field dependence 
of the field induced electronic specific heat A7(H). The non- 
linear field dependence is very clear. The red solid line is a 
fit to the 



(95% normal state resistivity p„), and the zero resistiv- 
ity was achieved at about 6.8 K. By applying a magnetic 
field, the resistive transition broadens. Taking different 
criterions of resistivity we determined the upper critical 
field Hc2(95%p„), Hc2(50%p„), and the irreversibility line 
Hirr(0.1%pn)- It is clcar that there is a large difference 
between the Hc2(95%/3„) and Hirr(T). We will argue that 
this may be induced by the strong SC fiuctuations. 

The raw data of specific heat was shown in Fig. 4(a). A 
SC anomaly appears at about 7.6 K. Since the Ru has a 
Tc at 0.49 K and a quite small normal state specific heat 
coefficient (7^" = 2.8mJ/molK^), a slight correction of 
about 0.42 mJ/molK^ was made to the data. By ap- 
plying a magnetic field, the SC anomaly shifts to lower 
temperatures. It is interesting to note that the transi- 
tion is not shifted parallel down to the low temperatures 
(the so-called mean- field like), rather the SC anomaly is 
suppressed. This kind of suppression was clearly seen in 



the cuprate superconductors Pro ««LaCen i9CuOa-a|14| 
and was ascribed to a strong SC fluctuation. Combin- 
ing with the resistive broadening under a magnetic field, 
we would argue that there is also a strong SC fluctu- 



ation in LaRu3Si2. As for a three dimensional system 
judged from our band structure calculations, this kind 
of strong SC fluctuation may suggest that the super- 
fluid density is low. Another interesting point shown in 
Fig. 4(a) is that the normal state specific heat (SH) coef- 
ficient C/T shows a non-Unear dependence on T^ down 
to the lowest temperature (0.38 K). This is clearly de- 
viating from the prediction of the Debye model. Tak- 
ing the slope of C/T vs. from the low tempera- 
ture data, we get the Debye temperature To ~ 284 K. 
The phonon contribution calculated based on the De- 
bye model Coebye oc {T /Tof j^"'^ [e^e^ / {e^ - \)'^]dx 
is shown by the red dashed line. One can see that the 
Debye model is seriously violated. It is naturally ques- 
tioned whether this violation is induced by some elec- 
tron correlation effect. For a non- Fermi liquid with three 
dimensionality, the enhanced electron-electron interac- 
tion will give an extra contribution to the electronic 
specific heat[ll C^-e = -AT'^lnT with n = 1 to 3. 
This gives the correction to the Fermi liquid descrip- 
tion, n = 1 corresponds to the case of strong correlation, 
like in Heavy fermion systems [l^, while n = 3 corre- 
sponds to a weak correlation. Thus we fit the data with 
C = 7„T + /3T3 - AT"-lnT and found a very good fitting 
when n takes 2 and 7„ = 36.i8mJ/molK^, as shown 
by the sohd fine, leading to (3 — lAl&mJ / molK^ and 
A = 3MmJ/molK^. Therefore we intend to conclude 
that the electron correlations may play an important role 
in the system. In Fig. 4(b), we derived the electronic spe- 
cific heat by subtracting the normal state background 
measured at 5 T. One can see that the low-T part of 
Ce/T exhibits a flat feature, indicating a full SC gap. 
There is a small upturn of Ce/T in the low-T region 
when the field is very weak, see the data of 0.25 T and 
0.5 T, which is attributed to the Schotkky anomaly due 
to the paramagnetic centers. Using an integral based on 
the BCS formula for electronic specific heat, we fit the 
data at ^qH — 0, 0.25 and 0.5 T and show with the sohd 
lines. In this way we obtained the data in the zero tem- 
perature limit for each field. For higher magnetic fields, 
it is known that the Schottky anomaly becomes weaker, 
we can determine the low-T data directly from the exper- 
iment data. Fig. 4(c) shows the temperature dependence 
of the entropy calculated using S — J^^ Ce/T'dT', it is 
clear that the entropy is conserved at Tc as judged by 
Ss — SnIt^ = 0, where Ss or Sn are the entropies of the 
superconducting state and the normal state integrated 
up to T, and Ss - Sn = J^iCs/T' - CN/T')dT'. After 
obtaining the low-T data, we derived the electronic SH at 
different magnetic fields and plot them in Fig. 4(d). Inter- 
estingly, a non-linear field dependence can be easily seen 
here. Further analysis finds that this non-linear depen- 
dence is actually different from the \/H relation (shown 
by the solid line in Fig. 4(d)) predicted for a clean super- 
conductor with line-node gap, the data below about 0.5 
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FIG. 5: (color online) (a) The energy bands obtained from 
the DFT calculations. The dense bands near Ef are derived 
from the Ru 4d orbitals. (b) The electronic density of states 
(DOS) from the band structure calculations. It is found that 
the DOS at the Ef are mainly contributed by the Ru orbitals. 
The DOS from the La and the Si atoms at Ef are negligible. 



T seems to be more linear. Our result here is certainly 
different from a linear relation as predicted for a single 
isotropic SC gap. A multigap feature would already ex- 
plain the data, but as seen from the flattening of Cg/T 
at T — > 0, we would argue that it is not the multigap 
feature, but the gap anisotropy that leads to the non- 
linear field dependence of Cg/T . A momentum resolved 
measurement is highly desired to uncover this puzzle. 

In order to have a comprehensive understanding to the 
data, we did the density-functional theory calculations by 
using the WIEN2k package jl6j utilizing the generalized 
gradient approximation [l7| for the exchange-correlation 
potential. As shown in Fig. 5(a), the bands around Fermi 
level are mainly contributed by Ru 4d. The Si 3p bands 
are very wide and have some hybridization with Ru 4d. 
Further analysis of the calculation shows that the crystal- 
field splitting upon Ru 4d orbitals is quite weak, con- 
sequently all Ru 4d electron should play an important 
role in the electron conduction and related supercon- 
ductivity. There are several bands crossing the Fermi 
level, which leads to complicated 3D Fermi surfaces, this 
will be presented elsewhere. Since the band closed to 
Fermi level is narrow, and the density of state at Fermi 
level is high as shown in Fig. 5 (b). We also perform 
spin polarized calculations to check the possible mag- 
netic instability. The calculation shows that the ferro- 
magnetism is not stable for this compound. We cannot 
find any strong nesting effect in the Fermi-surface, thus 
the SDW order is unlikely. Worthy of noting is that all 
the five Ru 4d orbitals contribute to the conduction in 
LaRu^Si2, which is very similar to the case of the iron 
in the iron-pnictide superconductors [isjl . Actually a Ru- 
based compound, namely LaRu2P2, is a superconductor 
with Tc ~ 4.1 K, which has the similar structure of the 
BaFe2~.xCoxAs2 superconductor (loi [20j . and probably 



they share the same superconducting mechanism. This 
reminds us that the correlation effect may play some roles 
in the superconductivity of LaRu^Si2- 

In summary, resistivity, magnetization and specific 
heat have been measured in a Ru-based superconduc- 
tor LaRu-iSi2 with Tc of about 7.8 K. The temperature 
dependence of the normal state specific heat coefficient 
C/T deviates clearly from the Debye model, but shows 
the possible evidence of electron correlations. The su- 
perconducting transitions measured by both resistivity 
and specific heat self-consistently present the evidence 
of strong superconducting fluctuations, resembling that 
in the cuprates. The field induced quasiparticle density 
of states show a non-linear magnetic field dependence, 
which is argued as a gap anisotropy. Combining the nov- 
elties found both in the normal state and the supercon- 
ducting state, we argue that the electron correlations may 
play an important role in the occurrence of superconduc- 
tivity in LaRu^Si2- 
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